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^vj ■ We present an analysis of the constraints on the amplitude of primordial non-Gaussianity of local 

type described by the dimensionless parameter /nl- These constraints are set by the auto-correlation 

oW, functions (ACFs) of two large scale structure probes, the radio sources from NRAO VLA Sky Survey 

^ . (NVSS) and the QSO catalogue of Sloan Digital Sky Survey Release Six (SDSS DR6 QSOs), as 

well as by their cross-correlation functions (CCFs) with the cosmic microwave background (CMB) 

temperature map (Integrated Sachs- Wolfe effect). Several systematic effects that may affect the 

observational estimates of the ACFs and of the CCFs are investigated and conservatively accounted 

for. Our approach exploits the large-scale scale-dependence of the non-Gaussian halo bias. The 

derived constraints on /nl coming from the NVSS CCF and from the QSO ACF and CCF are 

weaker than those previously obtained from the NVSS ACF, but still consistent with them. Finally, 

we obtain the constraints on /nl = 53 ± 25 (Icr) and /nl = 58 ± 24 (Icr) from NVSS data and 

r-j ' SDSS DR6 QSO data, respectively. 
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C/5 , I. INTRODUCTION 



The physical mechanisms responsible for the generation of primordial perturbations seeding present-day large-scale 
vN ' structure, may leave their imprint in the form of small deviations from a Gaussian distribution of the primordial 
perturbations. Searches for primordial non-Gaussianity can thereby provide key information on the origin and evo- 
lution of cosmological structures (e.g., ref. [l[ and references therein). Although the standard single-field, slow- roll, 
qC canonical kinetic energy and adiabatic vacuum state inflation generates very small non-Gaussianity, any inflationary 

^—{ ' model that deviates from this may entail a larger level of it (refs. 0,3 and references therein). 

J ' ' Deviations from Gaussian initial conditions are_often taken to be of the so-called local type and are parameterized 

^— N by the constant dimensionless parameter /nl 



* = + /nl (</-'- (</-'» , (1) 

where $ denotes Bardeen's gauge-invariant potential and (/) is a Gaussian random field. In the literature there are 

K> , two conventions: in the large scale structure (LSS) convention $ is linearly extrapolated to z = 0, while in the cosmic 

5^ microwave background (CMB) convention it is evaluated deep in the matter era. Thus, /^l^ — [5(2 = oo)/(7(z — 

F^ , 0)]/™^ "^ I-^/nl^^j where g{z) denotes the A-induced linear growth suppression factor. In this paper we will use 

the CMB convention. 

A new method y, 10] for constraining non-Gaussianity from large-scale structure surveys exploits the fact that the 
clustering of dark matter halos -where galaxies form- on large scales increases (decreases) for positive (negative) /nl- 
In particular, a non-Gaussianity described by eq. (P), introduces a scale-dependent boost of the halo power spectrum 
proportional to l/k'^ on large scales {k < 0.03/i/Mpc), which evolves roughly as (1 -K z). LSS surveys covering large 
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volumes are needed to access the scales where the signal arises (e.g., ref. [Tl| and references therein). Among the 
many currently-available tracers of the LSS, the radio sources from NRAO VLA Sky Survey (NVSS) [12] and the 
QSO catalogue of Sloan Digital Sky Survey Release Six (SDSS DR6 QSOs) [1^ are particularly interesting since they 
span large volumes extending out to substantial redshifts [1J|. Indeed these source samples were shown to provide 
tight constraints on primordial non-Gaussianity J15l4l7| . 

Extragalactic radio sources are uniquely well suited to probe clustering on the largest scales. Radio surveys are 
in fact unaffected by dust extinction which may introduce in the observed large-scale distribution spurious features 
reflecting the inhomogeneous extinction due to Galactic dust. Moreover, due to their strong cosmological evolution, 
radio sources are very rare locally, so that radio samples are free from the profusion of local objects that dominate 
optically selected galaxy samples and tend to swamp very large-scale structures at cosmological distances; thanks to 
the strong cosmological evolution, even relatively shallow radio surveys reach out to substantial redshifts. The NVSS 
[12| offers the most extensive sky coverage (82% of the sky to a completeness limit of about 3 mJy at 1.4 GHz) with 
sufficient statistics to allow an accurate determination of the angular correlation function, ^/{O) on scales of up to 
several degrees [l^ ITa. 

In a recent paper [I3] it was shown that the observed NVSS Auto-Correlation Function (ACF) hints at a positive 
value of /nl at about the 3 a confidence level. This is because the ACF is found to be still positive on angular scales 
9 > 4°, which, for the median source redshift (z„i ~ 1), correspond to linear scales where the correlation function 
should be negative if the density fluctuation field is Gaussian. A positive /nl adds power on large angular scales, 
accounting for the observed ACF. 

A cross-check for a positive /nl can be provided by an enhancement, compared to the Gaussian case, of the 
Cross-Correlation Function (CCF) of the CMB with LSS probes (late-time Integrated Sachs- Wolfe (ISW) [20J effect). 
Earlier analyses [H did not find any evidence for /nl > from the NVSS-CMB CCF, although the NVSS-CMB cross- 
spectrum showed an anomalous point, suggestive of some unaccounted systematic effect in the spatial distribution 
of NVSS sources and/or in the adopted CMB map. In this paper we will revisit this issue investigating possible 
systematics that may affect the results. 

We will also revisit the ACF of SDSS DR6 QSOs [H and their CCF. This was previously analyzed by ref. [21^ to 
exploit the high redshift regime for early dark energy models and for the ISW effect. Here we use up-to-date CMB 
maps to derive constraints on primordial non-Gaussianity, following the approach of ref. |15j . Optically selected 
QSOs are also well suited to test primordial non-Gaussianity, as they probe large-volumes and high redshifts and are 
not seriously affected by dust; their major contaminating source being stars from our own galaxy. Ref. p^] used 
an extension of the SDSS DR3 QSOs sample to constrain /nl- Here we revisit the analysis using an improved bias 
model, improved QSO catalog, updated knowledge of the sources redshift distribution, up-to-date CMB maps and 
complementary analysis methods. 

The structure of the paper is as follows. In §|n]we review the effects of primordial non-Gaussianity on the ACF and 
on the CCF, and the theory of the late-time ISW effect. Section Hill contains the analysis of ACF and CCF for NVSS 
radio sources and SDSS DR6 QSOs. In iJ IIVI we present the method used to derive constraints on /nl- Section IVl 
contains our main results. We conclude with a discussion and comparison with related work in iJ IVIl 

II. THEORETICAL FRAMEWORK 

In this section we briefly present the equations describing the ISW effect, and review the impact of primordial 
non-Gaussianity on cosmic observables relevant to our analysis, namely the halo mass function and the halo bias. 

A. Integrated Sachs- Wolfe effect 

The temperature anisotropy due to the ISW effect can be expressed as an integral of the time derivative of the 
gravitational potential $ over conformal time rj 

— (n) = -2 / $[?/, n(77o - il)]di] . (2) 

A CMB photon falling into a gravitational potential well gains energy, while loses energy when it climbs out of it. 
These effects exactly cancel out if the potential $ is time independent, such as in the matter dominated era, when the 
fractional density contrast 6„i is proportional to the scale factor a, so that $ = and no ISW is produced. However, 
when dark energy (or curvature) becomes important at later times, the potential evolves as the photon passes through 
it and <& ^ 0, producing additional CMB anisotropics. 



The late-time ISW effect can thereby be a powerful tool for probing dark energy and its evolution. However, the 
main contribution of the ISW effect to CMB anisotropics occurs on large scales that are strongly affected by cosmic 
variance. This problem can be overcome by cross-correlating the CMB temperature fluctuations with the distribution 
of extragalactic sources. 

The number density contrast in a given direction ni is: 



(5g(ni) 



f{z)dm{ni,z)dz 



/ bg{z) — {z)Sm{ni, z)dz , 



(3) 



where bg{z) is the scale-independent bias factor relating the density contrast of visible objects, (5g, to the mass density 
contrast ^,„ (^g — 6g(5m), and dN/dz is the normalized redshift distribution of the survey. 
The ISW temperature anisotropy in a direction n2 is: 
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The ACF of a complete sample of sources, C^^{0), and its CCF with a CMB map, C^'^{9), as a function of the angular 
separation 9 between ni and n2, can be written as: 
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where the P^ [cos(6')] are the Legendre polynomials, 9^ is the smoothing scale of the CMB map. Here Cf^ and Cf 
are the auto-correlation and cross-correlation power spectra, respectively, given by: 



Cj 



k^dkP{k)[If{k)] 



Cf = - / k'dkP{k)l]t'^{k)If{k) , 



in terms of the present day matter power spectrum, P{k), and of the functions If{k) and if'^ (k): 

If{k) = bg{z) — {z)Smik,z)je[kxiz)]dz, 
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je{x) being the spherical Bessel functions, and x the comoving distance. We use the publicly available package 
CAMB_sources^ 22] to calculate the theoretical ACF and CCF. 

B. Effects of Primordial Non-Gaussianity 

In the presence of primordial non-Gaussianity, the mass function n-t^Q{AI, z, f-t^-^) can be written in terms of the 
Gaussian one nQ^{M,z), for which a good fit to the results of simulations is provided by e.g., the Sheth-Tormen 
formula [2j], multiplied by a non- Gaussian correction factor |24l - [26| ^: 
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(11) 



where the normalized skewness of the density field, S's^m, is ex; /nl, ctm denotes the rms of the dark matter density field 
linearly extrapolated to z = and smoothed on scale R corresponding to the Lagrangian radius of a halo of mass M, 



^ Available at |http://camb. info/sources/ 1 

^ Although attempts have been made to derive directly an expression for the non-Gaussian mass function |27l . [28 
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FIG. 1: Effects of non-Gaussianity on the auto-correlation and cross-correlation power spectra (left panels), and on the ACF 
and CCF (right panels) for three different Gaussian halo bias models: b = 1.7 (bottom black curves), b = 1.1 4- 0.6/D{z) 
(middle red curves), fe = 1.1 -|- 0.6/D^{z) (top blue curves). Solid curves: /nl = 0; dashed curves: /nl = 100. 



and ^cc is the critical density for ellipsoidal collapse, calibrated on N-body simulations [29]. For high peaks (^oc/cm 3> 
1) and small /nl, <^oc is slightly smaller than the critical density for spherical collapse, Sc{z) — Ac{z)D{0)/D{z) where 
D{z) is the linear growth factor, and Ac(z) ~ 1.68 evolves very weakly with redshift. 
The large-scale halo bias is also modified by the presence of non-Gaussianity [3, llfl, l2a| : 



b-Moiz) - boiz) ~ 2{bG{z) - l)/NL^oc(2)aM(fc) 



(12) 



where the factor auik) encloses the scale and halo mass dependence. In practice, we find that, on large scales, 
auik) oc 1/fc^, is independent of the halo mass. 

We assume that the large-scale linear halo bias for the Gaussian case is |23| : 
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where z^ is the halo formation redshift, and Zq is the halo observation redshift. As we are interested in massive haloes, 
we expect that Zf ~ Zq. Here, q = 0.75 and p = 0.3 account for the non-spherical collapse and are a fit to numerical 
simulations (see also refs. [30 32]). 

Finally, the weighted effective halo bias is given by 



foNcl-^^min, Z, k, /nl) = 



Jm ■ rif^cdM 



(14) 



Afmin being the minimum halo mass hosting a source of the kind we are considering. 

Two things should be stressed when using eq. P^ : a degeneracy between be and /nl is expected at a given scale 
(the same amount of non-Gaussian bias can be given by different pairs of &g and /nl values; strictly speaking Bq is 



not a free parameter here, and the degeneracy is between Mmin, which is a free parameter, and /nl; however 6g is 
strongly dependent on M^i^); the 1/k'^ scale-dependence implies that large-scales are primarily affected by /nl, while 
small scales are mainly affected by be (see a quantitative discussion for this in ref. [13] )■ 

From the two left-hand panels of figure [2 one can see that a positive /nl enhances the amplitudes of the auto- 
correlation and cross-correlation power spectra especially on large angular scales {£ < 100). Consequently, the 
amplitudes of the ACF and CCF are also enhanced (right-hand panels of figure [T]). The ACF remains positive up 
to angular scales 6 > 4° where, for the adopted redshift distribution and a redshift-independent My^in, the ACF is 
expected to become negative. The lower left-hand panel of figure [T] shows that the effects of non-Gaussianity on the 
cross-correlation power spectrum become very large for the lowest multipoles. Actually, most of the non-Gaussian 
contribution to the CCF comes from £ < 10, corresponding to angular scales > 18°. Its detection therefore requires 
uniform sky coverage up to very large scales, but this is very difficult to achieve observationally. For example, in 
the case of NVSS sources, very small modulations of the source surface density on very large scales, such as those of 
interest here, are easily swamped by systematic variations of the survey effective depth correlated with the varying 
observing conditions [l3|. Also, the "integral constraint" [3j| needs to be dealt carefully. In fact, the ACF and CCF 
estimators involve differences between local values of some quantity and their mean value over the considered area 
(see eqs. (|16p and (fT7|) ). But if such area is not much larger than the correlation scale, the computed mean may differ 
from the true mean by a constant offset c, not known a priori, that must therefore be considered as a free parameter. 

III. OBSERVED ACF & CCF DATA 

In this section we describe two LSS probes, NVSS radio sources and high redshift SDSS DR6 QSOs, and present 
observational determinations of their ACF and CCF with the CMB. 

A. NVSS Radio Sources 

We start by confining our analysis to NVSS sources brighter than 10 mJy, since the surface density distribution of 
fainter sources suffers from declination-dependent fluctuations [Ig. Density gradients in the NVSS catalog become 
increasingly unimportant as the source brightness threshold is increased (see refs. jlSl. |34|). Also we mask the strip 
\b\ < 5°, where the catalog may be substantially affected by Galactic emissions. The NVSS source surface density 
at this threshold is 16.9 deg~ and the redshift distribution at this flux limit has been recently determined by ref. 
[35]. Their sample, complete to a flux density of 7.2 mJy, comprises 110 sources with Siaghz ^ lOmJy, of which 78 
(71%) have spectroscopic redshifts, 23 have redshift estimates via the K-z relation for radio sources, and 9 were not 
detected in the K band and therefore have a lower limit to z. We adopt here the smooth description of this redshift 
distribution given by ref. [36| : 

dN/dz = 1.29 -I- 32.37Z - 32.892^ + 11. 13^^ - 1.25z^ . (15) 

The mean redshift of this sample is (z) = 1.23. From the NVSS catalog we construct a pixelized map using the 
HEALPix software package [37] with resolution A^side = 64, yielding pixel areas of 0.92° x 0.92°. 
Our ACF estimator w{9) reads: 

^(,)^ i^Ii^^^MJ^i^M, (16) 

where rij is the number of sources in the i-th pixel, /, is the un-masked fraction of the pixel area, n is the expectation 
value for the number of objects in the pixel area [2l| . The sum runs over all pixel pairs whose centers are separated by 
an angle 9. The equal weighting used here is nearly-optimal because of the uniform NVSS sky coverage and because 
on large scales the noise is dominated by sample variance. We use A'b = 9 angular bins, 1° wide. The first one is 
centered at ^ = 40' (on smaller angular scales the dominant contribution to the angular correlation function comes 
from multiple components of individual radio galaxies); the other 8 are centered at 1° ■ • • 8° (on still larger scales the 
signal may be affected or even dominated by spurious density gradients [l8l.[l9j'). 

The estimated NVSS ACF was previously found to be positive over the full range of angular scales we consider [Vq 
(see also ref . 181 ). although the integral of 'w{9) over the full survey solid angle vanishes by construction. As pointed 
out by ref. [3^] , for a Gaussian distribution of primordial density fluctuations and a realistic redshift distribution for 
NVSS sources, the w{9) must vanish at 6* ~ 4° and become negative on larger scales unless the redshift dependence of 
the bias parameter for radio sources is drastically different from that of optical QSOs. Consistency with the clustering 



properties of the latter sources can be recovered allowing for a small amount of non-Gaussianity (/nl = 62 ± 27 (1 a) 

Since the NVSS covers most of the sky it can safely be assumed to provide a "fair sample" of the universe, so that 
the derived angular correlation function is not affected by the "integral constraint" significantly. Nevertheless, in rcf. 
|17| we have explored whether the indications of an excess (compared to the Gaussian case) positive contribution to 
the ACF on large scales may be due to a unexpectedly large positive offset c (a negative c would obviously exacerbate 
the discrepancy with expectations from the Gaussian case). To this end we have added to w{9) a constant c and have 
marginalized over it, allowing this quantity to vary in the range [10~^, 10""*], where the lower limit means essentially 
zero (we work on a logarithmic scale) while the upper limit corresponds to the case where the offset accounts for 
the full clustering signal on the largest scales (see ref. [3^). The best-fit value of c is about ~ 10~^, i.e. negligibly 
small, confirming that the data do not show indications that w{9) is appreciably affected by the integral constraint. 
On the other hand, c is weakly constrained by the data, and values much larger than the adopted upper limit are 
formally allowed. Such large values of c would imply that Wtruo(^) = Westimated(^) — c is negative on scales 6* > 4°, 
and the indications of non-Gaussianity would no longer be statistically significant. Such values however are not 
physically plausible because c is necessarily negligibly small when, as in the NVSS case, the survey area is so much 
larger than any other relevant scale (see Sec. 9.4.2 of ref. [33|), even in the non-Gaussian case. Allowing also for the 
correction proposed by ref. [39] to account for the infrared divergence of the non-Gaussian halo correlation function 
the constraints on /nl, become /nl = 58 ± 28 (1 tr) [171. 

A potentially trickier issue are large-scale surface density gradients due to instrumental effects that may spuriously 
enhance the ACF estimates. This issue has been extensively discussed by ref. [l^ (see also ref. [13]) who concluded 
that this effect is negligible if we restrict ourselves to a flux limit of 10 mJy, as we did. An upper limit to the 
magnitude of the effect is an offset of 10"'*, that, as noted above, would account for the estimated clustering signal 
on the largest scales. Having allowed c to take on this value, we have automatically allowed for the corresponding 
decrease of the statistical significance of /nl > 0. As a further test of the possibility that the large-scale clustering is 
spuriously enhanced by surface density gradients due to instrumental effects we have redone the analysis restricting 
ourselves to sources above 20 mJy, for which such effects are negligible on all scales, so that c = 0. In spite of the 
poorer statistics we find indications of /nl > at approximately the same significance level as for the 5 > 10 mJy 
sample: /nl = 73 ± 32 (1 cr error) and 19 < /nl < 139 (95% confidence level). 

To measure the CCF between the NVSS number density map and the CMB map, we use the following estimator: 



£'^iO) = ^Y.(Ti-T)'^^J^, (17) 

where Tj is the CMB temperature in the i-th pixel and T is the mean (monopole) value for the CMB temperature 
in the unmasked area. The results presented here rely on the 7 years Internal Linear Combination (ILC) map of the 
CMB provided by the Wilkinson Microwave Anisotropy Probe (WMAP) team [40]. We have checked that the using 
other CMB maps (the 5 years ILC map by the WMAP team and the improved (cleaner) map by ref. [4l|) does not 
change the results in any significant way. We adopt the WMAP KQ75 mask, excluding about 30% of the sky at 
low Galactic latitude, to avoid most of the residual Galactic contamination. This mask is then combined with the 
NVSS sky coverage; regions of the NVSS catalogue that may introduce spurious features because of missing snapshot 
observations and over-dense regions associated with bright of extended sources (Table 1 of ref. [T^) are also masked. 

As shown in iJ IIIBI the effect on the CCF of primordial non-Gaussianity shows up primarily on very large angular 
scales, at variance with the ACF for which the non-Gaussian signal is localized a.t 9 — 2° — 5°, as shown by ref. 
[l7| . Thus the residual large-scale systematic fluctuations in the NVSS source density with declination due to the 
varying projection of the beam and the change in observing configuration, that ref. [18] found to be still present at 
the 10 mJy flux limit, may seriously bias the CCF estimate. While it is true that these fluctuations should be, on 
average, uncorrelated with the CMB signal and therefore, on average, should not contribute to the CCF signal, we 
deal here with a single realization of the Universe, and few independent modes; in practice, the jackknife procedure 
we employ to estimate the errors is not suited to deal properly with such an effect. We therefore restrict ourselves 
to 5* 1 .4(- !Hz > 20mJy, where such fluctuations are negligible. The normalized redshift distribution for the sample of 
ref. [33 does not change significantly if the fiux limit is increased to 20 mJy, and therefore can still be described by 
eq. p^ . The ACF for this sample does not differ in any significant way from that of the S'i.4ghz > lOmJy on the 
scales of interest. 

The CCF is computed in iVf, = 16 angular bins, spaced by 1° with 0° < 6* < 15°. The covariance matrix of the data 



points is estimated using the jackknife re-sampling method 42|, dividing the unmasked area into M = 30 patches. 
From the 30 ACF and the CCF estimates obtained neglecting a patch in turn, we compute the diagonal (variance) 
and off-diagonal (covariance) elements of the covariance matrix. Our estimate of the NVSS-CMB CCF, shown in 



figure [2l is fully consistent with previous estimates using different approaches [ij, [3J, |43I44 
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FIG. 2: Estimated CCF between NVSS radio sources brighter than 20 mJy and the WMAP7 ILC map. Error bars are jackknife 
estimated. The black sohd hne is the best fit non-Gaussian model, while the blue dashed line refers to the Gaussian case. In the 
latter case, Aamp is a free parameter and the fit is obtained for ylamp ~ 1.14. It is indeed very close to the best-fit non-Gaussian 
curve, but there is no physical justification for ^amp 7^ 1. 



B. SDSS DR6 QSOs 

The SDSS DR6 QSO catalog released by ref. 13] contains about 10^ objects with photometric redshifts ranging 
from 0.065 to 6.075 over a total area of 8417 deg (~ 20% of the whole sky). We refer the reader to ref. [ij| for a 
detailed description of the object selection with the non-parametric Bayesian classification kernel density estimator 
(NBC-KDE) algorithm. We use the electronically-published table that contains only objects with the "good" flag 
with values in the range [0, 6]. The higher the value, the more probable for the object to be a real QSO (see § 4.2 of 
ref. [13] for details). Furthermore we restrict ourselves to the "uvxts=l", i.e. to QSOs clearly showing a UV excess 
which should be a signature of a QSO spectrum (in this case we have iVqso ~ 6 x 10^ QSOs). In order to minimize 
the effect of Galactic extinction on the observed QSO distribution, we mask regions with Ag > 0.18. 

We fit the redshift distribution dN/dz of the DR6 QSO sample with a function of the form: 



dN ._ 13 z 



7n-{- 



-exp 



(18) 



The best-fit values of the parameters are m — 2.00, /? = 2.20, zq ~ 1.62; the mean redshift of the sample is z ^ 1.49. 

In spite of the high efficiency of the selection algorithm adopted to define the SDSS DR6 QSO catalog, some 

contamination from UV-excess stars is unavoidable. Following ref. [46|, we model the observed ACF, c**, as the sum 



of contributions from QSOs, 



and from stars, c^^, plus an offset, e, arising from cross-terms: 
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(19) 



where a is the efficiency of the QSO catalog, i.e. the fraction of true QSOs. As shown by ref. [46|, the KDE 
classification technique is efficient enough for the offset e to be safely neglected. The contributions from stars and 
QSOs can be disentangled exploiting their different dependencies on 6. In fact refs. [46| and [2l| showed that 6^^ (9) 
keeps almost flat up to large angular scales, and is therefore expected to dominate the signal on scales of several 
degrees. We have checked that the contribution from contaminating stars to the QSO-CMB CCF can safely be 
neglected. 

As in the NVSS case we used the HEALPix software to pixelize the QSO map at the Ngidc = 64 resolution and the 
jackknife re-sampling method to estimate the covariance matrix among data points. Our estimates of the ACF for 
the SDSS DR6 QSO catalogue and of the QSO-CMB CCF, shown in figure [31 are consistent with previous analyses 
[2l|,|43|. Note that we do not split the SDSS QSO in different redshift bins using the objects' photometric redshift 
estimates, but we consider the projected ACF (and CCF) of the full sample. 
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FIG. 3: Estimated ACF of SDSS DR6 QSOs and CCF between their distribution and the WMAP7 ILC map. Error bars are 
jackknife estimated. We also show the predictions from our best fit model (red solid lines). In the left panel the black dashed 
line and the blue dash-dotted line show the estimated contributions from QSOs and from stars contaminating the sample to 
the global model ACF (red solid line). In the right panel we also show the the best fit model for the Gaussian case (blue dashed 
line) obtained letting Aamp to be a free parameter; the fit is obtained for ylamp = 1.86. It is indeed similar with the best-fit 
non- Gaussian curve, but there is no physical justification for j4amp 7^ 1. 



IV. METHOD & DATA ANALYSIS 

We perforin a global fitting of cosmological parameters, including /nl, for the data of Sec. IIIII including also the 
datasets described below, using the CosmoMC ■^ package [43], a Markov Chain Monte Carlo (MCMC) code, modified 
to calculate the theoretical ACF and CCF. We assume purely adiabatic initial conditions and a flat Universe, with 
no tensor contribution to primordial fluctuations. The following six cosmological parameters are allowed to vary with 
top-hat priors: the dark matter energy density parameter flch'^ G [0.01,0.99], the baryon energy density parameter 
ilhh^ £ [0.005,0.1], the primordial spectral index Ug £ [0.5, 1.5], the primordial amplitude log[10^''As] £ [2.7,4.0], the 
ratio (multiplied by 100) of the sound horizon at decoupling to the angular diameter distance to the last scattering 
surface 0s £ [0.5, 10], and the optical depth to reionization r £ [0.01,0.8]. The pivot scale is set at fcso = 0.05 Mpc~^ 
and do not consider massive neutrinos and dynamical dark energy. 

Besides these six basic cosmological parameters, we have three more parameters related to the ACF and CCF data: 
the non-Gaussianity parameter /nlj th^ minimal halo mass M,„in, two offset parameters c, one for the ACF and one 
for the CCF, allowing for the effect of the integral constraint. The ranges for /nl and Afmin are the same as in ref. 
[17| . Those for the offsets were chosen to be large enough to account for the large-scale signals. Specifically we have, 
for the NVSS: cacf £ [10-^ lO"*] (see ^MM, cccf £ [-0.5,0.5]; for the QSOs: cacf £ [-2 x 10-^,2 x lO'^], 
cccF £ [—0.5, 0.5]. All the best fit values of these parameters are well within these ranges and far from the boundary 
values (see ^ \V\ . 

In the QSO case we have, in addition, the efficiency of QSO classification a. Several authors also treated as an 
additional free parameter the ISW amplitude Aamp defined by (?'^{9) — Ai^^^pC^^ {6) , where c^""" and c^"^ are the 
observed and theoretical CCF. We note that values of j4ainp 7^ 1 have no physical meaning, and are indicative of some 
unrecognized problem in the interpretation of the data. Therefore we fix ^amp = 1 in all our calculations, except to 
estimate the significance of the ISW signal, in order to compare with previous works in which this parameter was 
used. 

The model ACF 08^(6*) and CCF c<^^{0) are compared with the observed values cSg(6') and CCF £^'^{9), respectively. 



Available at |http://cosmoIoglst.info/cosmomc/[ 
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FIG. 4: Marginalized one-dimensional posterior distribution of the non-Gaussianity parameter /nl from different NVSS data 
combinations: ACF alone (black continuous curve), CCF alone (red dashed curve) and ACF+CCF (blue dot-dashed curve). 



through the Gaussian likelihood function: 



£acf = {2Ti)-'''^[AetC,,r^'^ x exp 
/:ccF = (27r)-^/2j^g^^/_pi/2 ^ g^p 



y^ g^^(gf-cf)(gf-cf) 



«j 



^ c'^'{^-<f){^-^) 



ij 



(20) 
(21) 



where dj and C[j 



are the observed ACF and CCF covariance matrices. 

The following cosmological data are also included in the fit: i) power spectra of CMB temperature and polarization 
anisotropies; ii) baryonic acoustic oscillations (BAOs) in the galaxy power spectra; iii) SNIa distance moduli. 

To deal with the 7 years WMAP (WMAP7) CMB temperature and polarization power spectra we use the routines 
for computing the likelihood supplied by the WMAP team ,48:]. The WMAP7 data are exploited only to improve the 
constraints on the six basic cosmological parameters, not to constrain /nl- 

The BAOs [49| can, in principle, measure not only the angular diameter distance, Da{z)^ but also the expansion 
rate of the Universe, H{z). However, the limited accuracy of current data only allows us to determine the ratio 
between the distance scale defined by ref. |50| : 



A, 



{l + zfDl{z) 



-,1/3 



H{z) 



and the comoving sound horizon at the baryon-drag epoch, rJzd) (see ref. 
distance ratio rs(zii)/ Dy{z), rs{zii) have been obtained by ref. [49[: 



(22) 



|51j). Accurate determinations of the 



rsizd)/D,iz = 0.20) 
r,(zrf)/D,(z = 0.35) 



0.1905 ±0.0061, 
0.1097 ±0.0036. 



(23) 



We adopt these values as Gaussian priors. 

The SNIa data yield the luminosity distance as a function of redshift which provides strong constraints on the dark 
energy evolution. We use the Union compilation data (307 samples) from the Supernova Cosmology project [52| . 
which include the samples of SNIa from the (Supernovae Legacy Survey) SNLS and from the ESSENCE survey, and 
span the redshift range < z < 1. 55. In the calculation of the likelihood from SNIa we marginalize over the nuisance 
parameter as done in refs. 53l |54|. 

Furthermore, we add a prior on the Hubble constant, Hq = 74.2 ± 3.6 km/s/Mpc given by ref. [55|. Finally, in the 
analyses of NVSS sources and SDSS QSOs we set the minimal halo mass at Mmin > 10^^h~^MQ consistent with the 
results of ref. [563 . 
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FIG. 5: Cross-correlation between the ILC CMB and a Galactic emission template comprising synchrotron, free-free and dust 
contributions (black squares). The Galactic template is our best guess of the residual Galactic contamination in the ILC 
obtained as described in ref. 58\. Error bars are jackknife-estimated. 



V. NUMERICAL RESULTS 

A. NVSS Radio Sources 

To check the significance of the detection of the late-time ISW effect from the NVSS-CMB CCF, following ref. [lH, 
we compute the best fit value of the amplitude Aamp in the Gaussian case (/nl = 0) and the jackknife error. We 
obtain: 



A 



amp 



1.14 ±0.50 (1 (Terror) 



(24) 



so that the significance is ~ ^amp/^A: i-e. — 2.3 cr, broadly consistent with earlier analyses |43l - t45| and very well 
compatible with its physical value (unlike ref. [14|). 

The constraints on /nl, obtained setting Aamp = 1 and marginalizing over all the other parameters, are summarized 
in Table m The NVSS CCF turns out to be only weakly sensitive to the minimum halo mass. Therefore we fix it to 
the best fit value from NVSS ACF data alone, Mn,in = lO^^-^''h~^MQ. While the NVSS ACF yields a positive /nl at 
more than 2 a confidence level, /nl = 58 ± 28 (Icr), as previously found by ref. (o], the NVSS-CMB CCF provides 
much weaker constraints: /nl = 29 ± 48 (Icr). The best fit value for the offset is c ~ —0.05 for 5i.4ghz > 20mJy. 
To check the effect of residual systematic fluctuations in the NVSS source density, we have repeated the calculation 
for S'i.4GHz > lOmJy. We find that the amplitude of the CCF decreases, consistent with the cross-correlation being 
somewhat blurred (a similar effect can be gleaned also from figure 11 of ref. (34|), and the constraints on /nl become 
/NL = 6±37(la). 

Other systematic effects may affect the CCF estimates (for a parallel analysis on the ACF see ref. |17|). A first 
possibility is that the CMB map is contaminated by residual foreground emissions that increase the noise and thus 
swamp the signal we are looking for. The contamination from point sources below the detection limit is increasingly 
diluted for larger and larger pixel sizes. However we find that the CCF estimate does not change in any appreciable 
way when we use larger pixel sizes (A'sidc = 32 and A'sidc = 16). Note that increasing the pixel size we also decrease 
the relative contribution of CMB fluctuations generated at the recombination, which have much higher amplitude 
than those due to the ISW effect. 

The impact of residual contamination from diffuse Galactic emissions is assessed by cross-correlating the CMB 
map with the standard dust, free- free and synchrotron templates, used also by the WMAP team (e.g., ref. ^57;]), 
smoothed to 1 deg resolution. No indications of significant cross-correlations are found. However, residual foreground 
contamination will be due a certain mixture of the different Galactic emissions, thus it is not expected to strongly 
correlate with single foreground templates. We thus build a typical foreground residual map to be correlated with 
the CMB map, as described in the following. 

Ref. [58| shows how to predict the residual contamination after the application of a linear component separation 
technique. The method requires to have a model of the data and to know the weights of the linear mixture for 



11 



u 

.< -6 




20 40 60 80 100 120 140 




20 40 60 



100 120 140 



FIG. 6: Two-dimensional marginalized 1 & 2a contours in the (/nl, logioicACp]) and (/nl,cccf) from the NVSS ACF and 
CCF data. 



TABLE I: 1, 2(7 constraints on the primordial non-Gaussianity from different data combinations. We report the mean values 
and the Bayesian central credible interval, marginalized over all other parameters. 



Datasets 


Non-Gaussianity /nl 


WMAP7+BAO+SN 


la C. L. 


2cr C. L. 


NVSS Radio Sources 


+ACF 


58 ±28 


[16, 114] 


+CCF 


29 ±48 


[-50, 145] 


+ACF+CCF 


53 ±25 


[10, 106] 


SDSS DR6 QSOs 


+ACF 


42 ±28 


[-19,93] 


+CCF 


60 ±42 


[-20, 145] 


+ACF+CCF 


58 ±24 


[12, 94] 



component separation. 

Our data model is composed as follows: for the frequency dependencies of the components is based on model M2 of 
rcf. 59] , obtained from an analysis of the WMAP 3 years data; for the spatial morphology is based on the previously 
mentioned foreground templates. The component-separation linear-mixture weights are the ILC weights adopted to 
build the CMB map, available from the LAMBDA site. 

The foreground residual map we obtain is found to have a statistically significant negative CCF with the ILC map 
(figured]), indicating that the separation of the Galactic emission from the CMB map is not perfect. The amplitude 
of this spurious cross-correlation is however too small to affect significantly the NVSS-CMB CCF. 

If we combine with the NVSS ACF and CCF data we obtain the following constraint on /nl: 

/nl = 53±25 (Icr error), (25) 

or 10 < /nl < 106 at 95% confidence level. This result is compatible with previous estimates [iSl. l60l - [66| . and with 
the WMAP7 limits H. 

In figureini we plot the two-dimensional constraints on (/nl, cacf) and (/nl, cccf) from the NVSS ACF and CCF 
data. As noted above, a positive offset means that the true ACF or CCF is lower than the estimated one, thus 
weakening or spoiling indications for non-Gaussianity. On the contrary, a negative offset points to larger values of 
/nl- As discussed in section UlI Al the integral constraint should have a negligible effect on the NVSS ACF but an 
offset of up to 10^^ rnay be induced by large scale density gradients of instrumental origin. This value has been 
adopted as the upper boundary for the allowed range for cacf. 
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FIG. 7: Marginalized one-dimensional posteriors of the non-Gaussianity parameter /nl from different SDSS QSO data com- 
binations: ACF alone (black continuous curve), CCF alone (red dashed curve) and ACF+CCF (blue dot-dashed curve). See 
text for more details. 



B. SDSS DR6 QSOs 



As in Section FVAI to assess the significance of the detection of a cross-correlation between the distribution of SDSS 
DR6 QSOs and the CMB map we compute the best fit value of Aamp and its error keeping all the cosmological 
parameters fixed at the WMAP7 values. We find: 



^amp = 1.86 ±0.80 



(1 cr error) , 



(26) 



implying a ~ 2.3 cr detection of the late-time ISW effect, consistent with previous analyses. 

The constraints on /nl from the ACF, the CCF and their combination are given in Table U and illustrated in 
figure [71 As for NVSS sources, these constraints are obtained setting Aamp = 1 and marginalizing over all the other 
parameters. Allowing for the integral constraint may be very important for this data set, as the surveyed area is only 
~ 1/4 of the NVSS area and its effective angular size is not much larger than the scale over which positive correlations 
are induced by the /nl value indicated by the NVSS ACF data. As for the NVSS case we added constants c, as free 
parameters, to the ACF and to the CCF, and marginalized over them. In the ACF case, the effect of this constant 
is negligible. Its best fit value is cacf = —2 ■ lO"**, and we find /nl = 42 ± 28 (Icr); if we set cacf = we get 
/nl = 32 ± 19 (Icr). As expected, the impact of the constant offset is much more relevant for the CCF. We find 
ccCF = —0.15 and /nl = 60 ± 42 (1 a); setting cccf = we get /nl = ± 10 (1 cr). Clearly marginalizing over the 
cccF parameter has a big effect on the recovered /nl value. The larger the sky coverage of the survey the better this 
constant is determined and the less is its impact on /nl- 

Combining the ACF and CCF data we find /nl = 58 ± 24 (1 cr), cacf = -3 ■ lO"'', cccf = -0.13. The value of /nl 
is nicely consistent with that obtained from the analysis of NVSS data. The best fit value of the selection efficiency of 
the QSO catalogue, kept as a free parameter, is a = 97.1% ± 1%, consistent with the estimate by ref. [13] who claim 
an efficiency of over 97% for the UVX sub-sample considered here. 

Finally, in figure |8] we plot the two-dimensional constraints on (/nl, 10'* cacf) and (/nl, cccf) from the SDSS QSO 
ACF and CCF data. Both cacf and cccf are anti-correlated with /nl. A negative cacf and cccf enhances the 
estimate of /nl • 



VI. DISCUSSION & CONCLUSIONS 



All the datasets we considered here yield broadly consistent results. The NVSS ACF drives the positive signal 
16 < /nl < 114 at 95% confidence level. The SDSS QSO ACF and CCF individually yield results consistent with 
/nl = 0. This effect arises because of degeneracies in a high-dimensional parameter space. In particular, there is 
a degeneracy in the Mmin — /nl plane, with large values of Mmin corresponding to low /nl values, which, when 
one marginalizes over Mmin lowers the mean /nl value. The combination ACF+CCF however seems to break this 
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FIG. 8: Two-dimensional marginalized 1 & 2cr contours in the (/nl, IO^cacf) and (/nl,cccf) from the SDSS QSO ACF and 
CCF data. 



degeneracy (ACF depends on 6^ while CCF depends on b) cutting out the large Mmin tail and thus increasing the 
recovered /nl value. In addition the relatively large and negative cccf skews the maximum posterior /nl value for 
the CCF to larger positive values. Before drawing our conclusions we compare our findings with previous work in the 
literature where very similar data-sets were used. 



A. Comparison with Previous Works 

The work presented here is closely related to refs. [iajll3 ^^^ to the NVSS-CMB cross power spectrum as obtained 
by ref. [1J|. The most immediate difference in the analysis is that we work with the ACF and CCF while previous 
works has used the auto and cross power spectra. Correlation functions and power spectra are spherical harmonics 
transform pairs and should thus be fully equivalent. The two approaches however are effected by systematic effects in 
slightly different ways: whenever possible we try to directly compare the two approaches. In addition, compared with 
previous analyses, we marginalize over a suite of possible systematic effects, modeled by the constant offsets cacf and 
CCCF and over a possible stellar contamination and QSO efficiency parameter. 

Ref. [l^ used the NVSS-WMAP 5 year map cross power spectrum as provided by ref. [lj| to obtain /nl — 236±127 
(1 a). The main differences with our NVSS CCF analysis are: i) ref. [lg| used a constant bias, here we use eq. p^ . 
we include the non-Gaussianity correction in the mass function, thus have a free parameter given by the minimum 
host halo mass for the sources; ii) refs. |lj, [ig used WMAP 5 years map, while here we use up-to-date cleaned CMB 
maps; in) ref. [14] use NVSS sources much dimmer (> 2.5 mJy) than our sample (> 20 mJy). They have a higher 
source density but we have seen that density gradients are most important for dim sources. Finally, we do not attempt 
to model the redshift distribution of NVSS sources but use directly the one that has been measured; this results in 
our errors on /nl to be smaller, but the two results are fully consistent. 

Ref. [i3 obtained non-Gaussianity constraints using both the NVSS and SDSS QSO data. For their NVSS analysis 
they also used it only in cross correlation with the CMB (they do not use the NVSS sources power spectrum) and. 



like ref. [16|, used the cross-correlation data and errors of ref. [lj|. In their analysis the combinations: bias times 
source redshift distribution, b ■ dn/dz, and the bias as function of redshift, b{z), as well as /nl are determined from 
the correlation function data themselves, here instead we use the actual redshift distribution and use the minimum 
halo host mass as a parameter, rather than b{z = 0). In their reported ISW results several data sets are combined: 
the SDSS QSOs, the NVSS, 2MASS etc. although they report their results to be dominated by NVSS. Our NVSS 
CCF constraints on /nl are consistent with theirs although with smaller error-bars due to the fact that we do not 
need to estimate the source redshift distribution from the correlations themselves. 

For their SDSS QSO analysis ref. J5j considered both the auto and cross correlation signals. Ref. [151 used an 
extension of the DR3 QSO sample that include many of the sources that subsequently were released with the DR6 



sample: we use the more complete and better calibrated final official SDSS DR6 QSO catalog release [13|. Ref. [15 1 
found that their QSO sample at z < 1.45 seem to suffer from contamination which they described to systematic 
calibration errors, and discarded this sample from their analysis. We have checked for this effect and find no evidence 
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in the sample we use that the z < 1.45 QSOs have different contamination than the z > 1.45 or that the two samples 
give different estimates of /nl- Our sample is more than double the size than theirs. 

In our analysis we do not split the QSO in redshift bins according to each source photometric redshift estimate 
as ref. [15] do, but we consider the projected correlation of the full sample. Ref. [l^ had therefore to estimate the 
source redshift distribution from the (cross) correlation signals themselves for each of their sub-samples; we avoid this 
by using the independently-estimated full sample source redshift distribution. We also leave the stellar contamination 
as a free parameter to marginalize over, while ref. [15| kept it fixed at the fiducial value. 

The two analyses also use different bias models: ref. [l^ used a functional form b{z) ~ (1 -I- z)^ or b{z) ~ l/D{z) 
with different constant offset depending on assumptions about "recent merging" activity. We use the extended Press- 
Schecter bias formulation of eq. P^ . We have however verified that if we use their bias model, fix the QSO efficiency 
(i.e. fraction of stellar contamination in the catalog) to the fiducial value and do not marginalize over the constant 
offset we recover even central /nl values fully consistent with theirs (/nl = 6 ± 18 (Icr)). We therefore conclude 
that the largest effect driving the different central /nl values we recover is the treatment of systematic effects -stellar 
contamination and offset. 

We should stress here that for the data in common between the different analyses, the results we find are fully 
consistent (i.e., at the la confidence level) with both refs. [16] and jl^] results. Our positive /nl signal is driven by 
the NVSS auto-correlation function and by the combination of SDSS DR6 QSO ACF and CCF where the combination 
breaks the degeneracy with Afniin pushing he marginalized /nl values slightly up. 

Better knowledge of the bias of the sample and of systematic effects such as stellar contamination, calibration and 
selection effects are clearly needed to make further progress on this front; this is especially relevant to future data-sets 
where the large volume surveyed will further reduce the statistical errors. 

B. Conclusions 

We have investigated the constraints on the parameter /nl, characterizing primordial non-Gaussianity of the local 
type, coming from the currently available data sets best suited for this purpose, namely the observationally determined 
auto-correlation functions (ACFs) of NVSS radio sources and the SDSS-DR6 QSOs, and the cross-correlation functions 
(CCFs) of the surface density distributions of these sources with the WMAP 7 years CMB map. The mean redshift 
of NVSS sources is z ~ 1.23, and that of SDSS-DR6 QSOs is z ~ 1.49. 

As shown by ref. [l3l, a positive /nl value (/nl = 58 ± 28 (1 a), after allowing for the possible effect of the integral 
constraint), consistent with the constraints set by previous analyses, and especially by WMAP data [4^ (/nl = 32±21 
(Icr)), can explain the positive ACF signal detected on angular scales > 4° (see refs. [la [3), where the ACF is 
expected to be < for Gaussian primordial perturbations. The QSO ACF yields constraints (/nl = 42 ± 28 (1 a)) 
very similar to those inferred from WMAP 7 year data. 

A signature of /nl > should also be present in the CCF between tracers of large scale structure and the map 
of CMB anisotropics. Most of the signal, however, comes from low multipoles {i < 10) corresponding to very large 
angular scales where the signal may be blurred by systematic effects. In the case of NVSS sources we have indeed 
found that the amplitude of the CCF with the CMB increases somewhat as we increase the adopted flux limit from 
5'i.4GHz = 10 mJy, which is found to be adequate for the ACF analysis but is still affected by small spurious large scale 
surface density gradients, to S'i.4ghz = 20mJy where such gradients become negligible. Other possible systematic 
effects, such as contamination of the CMB map by point sources below the detection limit or by residual Galactic 
emissions have been shown to be of minor importance. The NVSS-CMB CCF however turned out to provide rather 
loose constraints on /nl (/nl = 29 ± 48 (1 a)). 

The limited area covered by the SDSS-DR6 QSO catalog implies that the CCF is liable to a quite substantial effect 
of the integral constraint, at variance with the NVSS case, where the area is 4 times larger (the ACFs of both NVSS 
sources and SDSS QSOs are found to be insensitive to this effect). We have allowed for the integral constraint by 
adding to the estimated CCF a constant offset, treated as a free parameter, and marginalizing over it. In this way we 
obtained /nl = 60 ± 42 (1 a). 

Finally, we obtain the constraints on /nl = 53 ± 25 (1 a) and /nl = 58 ± 24 (1 a) from NVSS data and QSO data, 
respectively. The tantalizing hints of a positive /nl reported by ref. [171] have survived the present more extensive 
analysis, involving other data sets. In order to reconcile the data with /nl = 0, we need to show that either the NVSS 
catalog is affected by many-times larger systematic effects (spurious large-scale correlations introduced by density 
gradients) than known and quantified so far, or that the redshift evolution of the bias of these sources is radically 
different from that of their radio-quiet, optical counterpart. 

Alternatively, one may question the accuracy of jackknife error-estimation especially when there are many non-zero 
off-diagonal elements in the covariance matrix. From ref. [67| . we infer that jackknife can underestimate parameter 
errors by up to 30%. If the error bars were to be increased, /nl would become more compatible with zero. Much 
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tighter and robust constraints on /nl should be provided by the forthcoming surveys with ASKAP [68] and MeerKAT 
^^, the SKA pathfinders, and by much larger area QSO surveys. 
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